Numerical Study of the Vortex Phase Diagram Using the Bose Model 

in STLS approximation 
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We study the phase diagram of the flux hues using the mapping to 2D bosons in the self-consistent- 
field approximation of Singwi, Tosi, Land, and Sjolander (STLS). The pair correlation function, 
static structure factor, interaction energy, and spectrum of the excited energies are calculated over 
a wide range of the parameters in this approximation. These quantities are used for studying the 
melting transition from the Abrikosov lattice into the entangled vortex liquid. The resulting B — T 
phase diagram is in good agreement with the known estimates for the vortex lattice melting and 
the Monte Carlo simulations. We also discuss the effect of van der Waals interaction, induced by 
thermal fluctuations, together with the repulsion potential on the phase diagram. 



I. INTRODUCTION 

The phases of flux lines (FLs) in high temperature su- 
perconductors are the subject of many current experi- 
mental and theoretical investigations [0j2j. In the classi- 
cal (Abrikosov) picture of type-II superconductors, FLs 
penetrate a clean material for fields H > Hci , to form a 
triangular solid lattice. This mean field phase diagram 
is modified by inclusion of the strong thermal fluctua- 
tions and play an important role in high-Tc materials 
|^|-^. The presence of various forms of disorder (such 
as point or columnar disorder), also affects on the mean 
field phase diagram. In this paper we typically focus on 
the phase diagram of the pure materials, however, there 
are many interesting experimental and theoretical works 
describing the effect of both thermal fluctuations and dis- 
order on the phase diagram [7|-pX|]. 

It is well known that the thermal fluctuations lead to 
melting of the vortex lattice and appearance of a vortex 
liquid. It has now become possible to observe the melt- 
ing transition experimentally. The indirect observations 
based on the resistivity measurements, and re- 
cent experiments p"3|-p^ based on the measuring a jump 
in the magnetization or the latent heat show a first-order 
vortex lattice melting transition. On the other hand, vor- 
tex lattice melting has been studied theoretically using 
various approximations. Early works using the renormal- 
ization group or density functional theory have 
indicated a first-order transition. Elastic theory com- 
bined with the Lindemann criterion produce a melting 
line in good agreement with the experimental observa- 
tions ]T9[ |. There are also a large interest in numerical 
simulations for studying this problem. An interesting 
work in this direction is done by Nordborg and Blat- 
ter 1 20 1 which present an extensive numerical study of 
the vortex matter using the mapping to 2D bosons and 
path-integral Monte Carlo simulations. 

It was suggested by Nelson that the vortex system is 
equivalent to a system of interacting bosons in two di- 
mensions (Bose model) This mapping predicts a 
melting transition into an entangled vortex liquid. There- 



fore, the problem of a vortex system maps to a system 
of N bosons in two dimensions interacting through the 
potential V{r) — g^Ko{r/\), where Kq is the modified 
Bessel function, A is the London penetration depth, and 

is a constant that scales the energy of interaction. In 
the language of the vortices, this potential comes from 
the interaction between vortices in the London theory, 
and g is related to the elastic moduli of the vortex lat- 
tice. The Bose model differs from the real vortex system 
(see next section) , however it still contains the main part 
of the interaction ||2^. Hence, it would be reasonable and 
interesting to use this model for describing the properties 
of the vortex phase diagram. More recently , the Bose 
model has been used in a numerical study of the vortex 
matter. Some physical quantities such as structure fac- 
tor and superfluid density in different temperatures are 
given addressed, and the first order vortex lattice melt- 
ing transition into an entangled vortex liquid is approved 
by numerical simulations. In the language of the boson 
system, this transition is related to the quantum phase 
transition from a Wigner crystal to a superfiuid. 

The Bose model idea also allows for using the many 
body techniques. In this work, we apply the self- 
consistent-field approximation of Singwi, Tosi, Land, and 
Sjolander (STLS) |^l[], and calculate the static structure 
factor, pair correlation function, interaction energy and 
the spectrum of the excited energies for different mag- 
netic field strengths and temperatures. The STLS ap- 
proximation has originally been proposed for describing a 
degenerate electron gas and has been used successfully to 
study a variety of other systems too [p2|j2^ ] . In the STLS 
theory, the correlation effects are incorporated through a 
static local-field correction, which is obtained numerically 
in a self-consistent way. We find numerical results for the 
static structure factor S{q) over a wide range of the pa- 
rameters. From the calculation of S{q), we present the 
results for the pair correlation function, the interaction 
energy and spectrum of the excited energies. Different 
behaviors of these quantities may be used for studying 
the phases of vortex lines. 

It is well known that the oscillatory behavior in g(r) 
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is a signature of the sohd phase. Therefore, the phase 
transition can be detected by looking at the behavior 
of the g{r) in a fixed temperature (magnetic field) but 
varying magnetic field (temperature). The solid-liquid 
transition can also be observed using the static structure 
factor. Disappearing of the peaks in the structure factor 
resembles the onset of the phase transition. Hence, us- 
ing the behaviors of the pair correlation function and 
static structure factor, the phase diagram can be ex- 
plained qualitatively, however, it is not possible to de- 
termine precisely the melting transition temperature. 

One of the transition signatures is the appearance of a 
special q on which the spectrum of the excited energies 
vanishes. So we have numerically investigated the disper- 
sion relation of the excited energies as a good quantity 
for the estimation of the B — T phase diagram and the 
melting temperature. Quantitatively our results for the 
excited energies are compatible with the expected results 
of the phase diagram and recent Monte Carlo simula- 
tions @. 

On the other hand, the direction dependent interaction 
for real vortices has very interesting consequences and 
predicts a van der Waals interaction even for straight 
vortices [Q. It is shown ||2j.|2^ that in the decoupled 
limit, 7 — !■ (where 7 is the anisotropy parameter), the 
van der Waals attraction is proportional to l/R^. This 
attractive interaction with entropic repulsion has very 
important outcomes for the low-field phase diagram of 
the anisotropic superconductors [^,^ . Recently Volmcr 
and Schwartz introduced a new variational approach 
to consider the effects of van der Waals attraction and 
repulsive interaction on the field phase diagram. We also 
consider the same model in the STLS approximation for 
studying the combined effects of the repulsive and attrac- 
tive potentials on the solid and liquid phases of the pure 
anisotropic or layered superconductors. 

The rest of this paper is organized as follows. In Sec. II, 
we shortly review the Bose model and discuss its appli- 
cability to the vortex system. The STLS approximation 
is briefly discussed in Sec. III. The numerical results for 
the repulsion interaction is presented and discussed in 
the Sec. IV. The results for the van der Waals interac- 
tion and its consequences in the phase diagram is given 
in Sec. V, and finally the conclusions appear in Sec. VI. 



where Ri (r) is a two dimensional vector representing the 
positions of the bosons, and T is the temperature of the 
system. In the Schrodinger picture the above action is 
equivalent to the two-dimensional Schrodinger equation, 
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(2) 



where the potential V{Rij) is proportional to the modi- 
fied Bessel function. 

It was pointed out by Nelson that the above ac- 
tion can be also interpreted as the London free energy 
for a system of vortex lines. The London free energy for 
a system of interacting vortices for a sample of length Lz 
is given by. 
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where e; « ^i'^eoao/ {2y/TTC), £0 = (*o/47^A)^ $0 is the 
quantum flux, ^ is coherence length, and oq is the lattice 
spacing. Comparing this functional free energy, which is 
referred as Bose model for the vortex system, with the 
action (^) shows the relationship between the parameters 
of the 2D bosons and vortices. 

The modified Bessel function Kq{R/\) describes a 
screened logarithmic interaction, Ko{x) ^ —ln{x) for 
X — *■ 0, and ^^^0(2;) ~ x~^/'^e~^ for a; ^ 00. Thus, the 
London penetration depth defines the interaction range. 
Note that according to the two-fluid model [^|J|] , the pen- 
etration depth diverges at zero-field transition tempera- 
ture Tc, and therefore the interaction range becomes con- 
siderably longer upon approaching Tc- 

As it is discussed in details in the Ref. ||2^ , in spite of 
the fact that the Bose model differs from the real vortex 
system in the choice of boundary conditions, the linear- 
alization leading to the elastic term, and the retarded 
interaction, contain the main parts of the interaction be- 
tween vortices and one expects that the results be in a 
rough quantitative agreement with those of real systems. 



III. STLS APPROXIMATION 



II. BOSE MODEL 



In the Feynman path- integral picture ]2^ , the system 
of N interacting bosons in two dimensions is described 
by the action 
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In this section we shortly review the principal equa- 
tions of the STLS approximation Q . These equations 
show the relation between the response function, static 
structure factor, and the local field correction. 

In the STLS approximation, the response function can 
be expressed as 



(4) 



In this equation xo is the response function of the free 
Bose gas and ip is the effective potential given by ■il;{q) = 
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v{q){l — G{q)), where v{q) is the Fourier transformation 
of the bare potential, and G{q) is local field correction. 
The STLS local field correction is 



nj (27r)^ q^ v{q) 



(5) 



where n is the density, and S{q) is the static structure 
factor and is related to response function as (?i = 1), 



S{q) 



(6) 



For a noninteracting two dimension Bose gas the free re- 
sponse function is given by. 
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where e(g) = q'^ / (2m) is the free particle energy, and is 
a positive infinitesimal quantity. 

Usiiig the Eqn. (^, one can calculate the integral 
Eqn. (^ analytically which leads to the following result 
for the static structure factor 
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The Eqs. (||) and (||) should be solved numerically for 
obtaining the S{q) self consistently. From the knowledge 
of the structure factor, the pair correlation function g{r) 
is calculated as 
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The interaction energy is also related to the structure 
factor as, 



values of the two parameters m and r^. For the boson 
system m can be considered as the mass of particles (in 
the unit of h = 1), and for vortex lines m is related to 
the temperature as m = SqEiX^ /T'^. Substituting the nu- 
merical values for parameters as Eq — 50K/A, 7 — 100, 
A = 1000 A, the equivalent temperature will be fixed as 
T w 500 / y/rri{K) . On the other hand, is the inverse 
of density corresponding to the particle density for boson 
system, or the density of FLs in the vortex matter. The 
mentioned relationship is expressed as i? = (/)o / (vrr^ ) , 
or B « 0.06/r2 (Tesla). 



One of our numerical results which was developed by 
iterating the Eqs. (||) and (^) is the behavior of the static 
structure factor which is shown in Fig. 1 for a fixed 
Ts = 0.6 {B 0.17(Tesla)) and different m's (temper- 
atures). 
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where vx{q) = Xv{q) and Sx{q) is its related static struc- 
ture factor. 

One can also compute the excited energy using the 
poles of Eqn. (0), leading to 



Loiq) = ^e^q) + 2ne{q)viq){l-G{q)). 



(11) 



In the next section, we present numerical results for 
various quantities of interest. 



IV. NUMERICAL RESULTS 



Fig. 1: The static structure factor in terms of ( 
and different m's. 



for 7-3 — 0.6 



Fig. 1 manifests the very expected behavior of the 
static structure factor, that is the peak of the S{q) is 
decreased by increasing T (decreasing m). Disappearing 
of the peak in the S{q) by increasing the temperature 
shows that the system undergoes a phase transition from 
the solid phase to the liquid phase. 



In this section we present the results of the numerical 
calculations for the interesting physical quantities. We 
numerically solve the set of Eqs. (o) and (g) with the re- 
pulsion potential defined in Eqn. (0), and find the static 
structure factor. The calculations are done for different 



We can also see the melting transition by fixing the 
temperature and increasing the magnetic field. Choos- 
ing m = 100 (T = 50K), we have found the S{q) for 
different values of r^, see Fig. 2. 



3 




rs: 


=0.4 


rs: 


=0.6 


rs: 


=0.8 


rs: 


=1 



10 



15 



the system is going to have a transition from the sohd 
phase to the hquid phase by increasing the temperature. 

We have also determined the interaction energy of the 
system using the Eqn. (|l0|). The results are shown in the 
Table 1. 



rs 


B 


Eint 


0.1 


6 


-2.15 


0.2 


1.5 


-1.78 


0.4 


0.38 


-1.36 


0.6 


0.17 


-1.09 


0.8 


0.09 


-0.89 


1 


0.06 


-0.74 



Table 1: The interaction energy, in STLS approximation for 
m = 50 and different densities (The unit of magnetic fields is 
Tesla). 



Fig. 2: The static structure factor for m = 100 and with 
different values of r^. 

The amplitude of the peaks becomes smaller with de- 
creasing the rs, while it disappears for high magnetic 
fields. So we will have the same scenario for describing 
the phase transition. 

The melting transition of the vortex lattice can also be 
discussed by studying the behavior of the pair correlation 
function. It is seen that the results of S{q) can be used 
for calculating the g{r) defined in the Eqn. (^. The re- 
suhs of g{r) for = 0.4 {B « 0.38(Tesla)) and different 
masses (temperatures) are plotted in Fig. 3. 



We observe that the numerical results show that the 
interaction energy increases by decreasing the strength 
of magnetic field. 

The results of the pair correlation function and static 
structure factor are in good qualitative agreement with 
the expected results of the phase diagram of the vortex 
system. However, for obtaining some quantitative de- 
scription, we use the spectrum of the excited energy of 
the system. We have plotted the uj{q) in terms of q, see 
Fig. 4. Fixing the magnetic field but for varying tem- 
peratures, the following graph resembles that io{q) tends 
to a zero value for a finite q . The appearance of mini- 
mum in LLj{q) is in correspondence with the fact that G{q) 
becomes greater than one for some ranges of m. 




Fig. 3: The pair correlation function for = 0.4 with dif- 
ferent temperatures. 

Fig. 3 shows that the pair correlation function has os- 
cillatory behavior for low temperatures, however, its am- 
plitude becomes shorter by increasing the temperature, 
and disappears for very high temperatures. The oscilla- 
tory behavior is a signature of the solid phase, therefore. 



Fig. 4: The dispersion relation of the excited energy for 
Ts = 0.6 for different values of m. 

The zero value of the dispersion relation indicate that 
the system undergoes a phase transition at that point. 
Because of numerical restrictions we are not able to find 
the exact value of the point where w(g) becomes zero. 
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However, we have realized by our numerical experience 
that when we are close to the transition point the q value 
of the minimum w{q) and the local field correction G[q) 
are not highly sensitive to the variations of m. There- 
fore, at least for determining the transition temperature 
we have ignored their dependence to m in the Eqn. (11). 
So we fixed the magnetic field and found the tempera- 
ture that becomes zero. The resulting B — T phase 
diagram is plotted in the Fig. 5. 



strength of the thermal vdW attraction. The function 
(j}{x) smoothly cuts off the power law part for R < X and 
is defined as, 



0, 



sin ^TT 



x-{xi+X2)/2 



X < Xi 

, Xi < X < X2 
X > X2 

(13) 




Fig. 5: The B ~ T phase diagram of the flux lines. The 
solid line is the predicted result in Ref. produced by Lin- 
demann criterion. Triangular points are the results of our 
numerical calculations. 

Proceeding to fit our data with the known result B = 
4cl(j)oeoei/{VST^) we found that the best fits 

can be achieved by fixing ~ 0.15 and it is observed 
that the results are supporting the reports of Monte Carlo 
simulations iHl. 



V. VAN DER WAALS INTERACTION 



with xi — 1 and X2 = 5. The choice of the cutoff function 
and the values of Xi and X2 is to some extent arbitrary 
1^ , ^ . The amplitude of the vdW attraction is given by 
dvdw ~ T/{2eodln^{Tr\/d)), where d is the layer spacing. 
Using the BSCCO numerical values for the parameters, 
one finds that Uydw ~ 2 x 10~^T/K which is 2 x 10~^ for 
the temperature T = IQOK close to critical temperature 



The potential defined in Eqn. (12) is applied for cal- 
culating the static structure factor in the STLS approx- 
imation. The results are plotted in the Fig. 6, where 

Avw = aydwVQ. 




TAking into account the direction dependent interac- 
tion for real vortices give rise to interesting results and 
predicts a van der Waals (vdW) interaction even for 
straight vortices [Q. Therefore, it would be interesting 
to consider the superposition of the short range attrac- 
tive and the long range repulsive interactions and study 
its consequence for the low-field phase diagram of the 
anisotropic superconductors. It is shown [E4 25 1 that in 
the decoupled limit, £ — > 0, the interaction potential is 
given by 

V{R) = vo (^Ko{R/X) - a„,„,/,(i?/A)^^ , (12) 

where, vq = 2eo/T measures the amplitude of the direct 
interaction between flux lines, and a^dw determines the 



Fig. 6: The static structure factor for fixed magnetic field 
rs = 0.2 with different temperatures and strengths of the vdW 
attraction. 

The results show that the importance of the vdW in- 
teraction with high magnetic field is valid only for the 
temperature close to the transition. 

One should note that the vdW interaction has very 
important consequences for low magnetic fields and tem- 
peratures. However, in our approach it is not possible 
to consider this case, and one has to use some other ap- 
proximations or Monte Carlo simulations for seeing these 
effects. 
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VI. CONCLUSIONS 

In this paper we studied the Bose model in STLS ap- 
proximation. We found the static structure factor, pair 
correlation function, interaction energy and spectrum of 
the excited energies for different values of the mass (tem- 
perature) and the density (magnetic field) . We discussed 
that how the results may be applied for the vortex mat- 
ter phase diagram. If one fixes the magnetic field (tem- 
perature) and calculate S{q) for different temperatures 
(magnetic field), the gradual disappearance of the peaks 
in S{q) manifests the existence of the phase transition. 
One can also use the pair correlation function for describ- 
ing the melting transition. We showed that the changing 
behavior of g{r) from oscillatory to a rather smooth one 
can be explored within our numerical scheme so that the 
hallmark of the transition from solid phase to the liquid 
phase can be observed . The results were in good quali- 
tative agreement with the phase diagram of the FLs. For 
estimation of the B — T phase diagram quantitatively, we 
invoke to the behavior of the excited energy spectrum, 
from which the resulting phase diagram supports the ex- 
pected results of the high temperature superconductors 
and Monte Carlo simulation quite well. 

We also added the van der Waals attractive potential 
and studied the effect of both repulsive and attractive 
potentials in the phase diagram. The results indicate 
that the vdW interaction for high magnetic field is only 
important for the temperatures close to the melting tem- 
perature. We emphasis that our approach doesn't work 
for the low magnetic fields and hence, it would be in- 
teresting to use some other methods and determine the 
effects of both repulsive and attractive potentials in the 
low magnetic fields and temperatures. 

To our knowledge this is the first time that the STLS 
approximation is applied for studying the vortex system, 
and our work show that STLS approximation is appli- 
cable for studying other aspects of the vortex systems 
and it might help to reveal some unknown properties in 
further investigations. 
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gari, J. Davoudi and M. Kardar. M. Kohandel acknowl- 
edges support from the Institute for Advanced Studies 
in Basic Sciences, Zanjan, Iran. We also acknowledge 
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and Mathematics, Tehran, Iran. 



[1] 
[2] 
[3] 

[4] 

[5] 

[6: 

[7 

[s: 

[9 

[lo: 
[11 

[12 

[13 
[14 
[15 
[16 

[17 

[18 
[19 

[20 
[21 
[22; 
[23; 
[24; 
[25 
[26; 
[27 
[28; 



G. Blatter et. al.. Rev. Mod. Phys. 66, 1125 (1994). 
E. H. Brandt, Rep. Prog. Phys. 58, 1465 (1995). 

A. I. Larkin and Y. N. Ovchinokov, J. Low Temp. Phys. 
34, 409 (1979). 

D. R. Nelson, Phys. Rev. Lett. 60, 1973 (1988); D. R. 
Nelson and H. S. Seung, Phys. Rev. B 39, 9153 (1989). 
M. P. A. Fisher, Phys. Rev. Lett. 62, 1415 (1989); D. S. 
Fisher, M. P. A. Fisher and D. A. Huse, Phys. Rev. B 
43, 130 (1991). 

M. Kohandel and M. Kardar, Phys. Rev. B 59, 9637 
(1999). 

D. R. Nelson and V. M. Vinokur, Phys. Rev. Lett. 68, 

2398 (1992); Phys. Rev. B 48, 1306 (1993). 

U. C. Tauber and D. R. Nelson, |cond-mat/9505002[ l, 
preprint (1995). 

D. Ertas and M. Kardar, Phys. Rev. Lett. 69, 929 (1992); 
ibid. 73, 1703 (1994); Phys. Rev. B 53, 3520 (1996). 
M. Kohandel and M. Kardar, preprint (1999). 

H. Safar et. al., Phys. Rev. Lett. 69, 824 (1992). 

M. Charalambous, J. Chaussy and P. Lejay, Phys. Rev. 
B45, 5091 (1992). 

E. Zeldov et. al., Nature375, 373 (1995). 

U. Welp et. al., Phys. Rev. Lett. 76, 4809 (1996). 

A. Schilling et. al.. Nature 382, 791 (1996). 

M. Roulin, A. Junod and E. Walker, Science 273, 1210 

(1996). 

E. Brezin, D. R. Nelson and A. Thiaville, Phys. Rev. B 
31, 7124 (1985). 

S. Senguptaet. al., Phys. Rev. Lett. 67, 3444 (1991). 
A. Houghton, R. A. Pelcovits and A. Subdo, Phys. Rev. 
B 40, 6763 (1989); G. Blatter, V. B. Geshkenbein, A. I. 
Larkin, H. Nordborg, Phys. Rev. B 54, 72 (199 6). 
H. Nor dborg and G. Blatter , 2ond-mat/9612025, preprint 



;ond-mat/9803041, preprint (1998) 



(1997); 

K. S. Singwi, M. P. Tosi, R. H. Land and A. Sjolander, 
Phys. Rev. 176, 589 (1968). 

K. S. Singwi and M. P. Tosi, Solid State Phys. 36, 177 
(1981). 

R. K. Moudgil, P. K. Ahluwalia, K. Tankeshwar and K. 
N. Pathak, Phys. Rev. B 55, 544 (1997). 
G. Blatter and V. B. Geshkenbein, Phys. Rev. Lett. 77, 
4958 (1996). 

A. Volmer, S. Mukherji and T. Nathermann, cond 



mat/9802090i preprint (1998) 



;ond-mat/980629g , preprint 



A. Volmer and M. Schwartz, 
(1998). 

R. P. Feynman, Statistical Mechanics, Addison- Wesely, 
Redwood City (1972). 

G. Blatter et. al., Phys. Rev. B 54, 72 (1996). 



2 



6 



